U. S. Department of Commerce 
National Bureau of Standards 



Research Paper RP1984 
Volume 42, May 1949 



Part of the Journal of Research of the National Bureau of Standards 

Transient Vibration in an Airplane Wing Obtained by 

Several Methods 1 

By Walter Ramberg 



Analysis of the flexural transients in a model airplane wing following a "soft" unsym- 
metrical two-point landing impact indicates that Williams' method, separating ' 'static" and 
vibrational response, is superior to the normal modes method used by Biot and Bisplinghoff, 
in leading to values in good agreement with experiment without using more than four modes. 
A method of coupled modes suggested by Levy leads to values of the response about as accu- 
rate as those from Williams' method with as few as three vibrational modes. 

None of the three methods showed rapid convergence for the response following a 
"hard" impact of the same shape as the soft impact, but with only one-fifth of the duration. 
Attempts to develop a traveling wave method for the solution of this problem have been 
unsuccessful so far. 



I. Statement of Problem 

The problem of determining the transients in 
an airplane following the landing impact has be- 
come important during recent years with the 
advent of large transport airplanes. 

It became apparent, as a result of several fail- 
ures in service, that the stresses developed during 
the landing of large transport airplanes could be 
far in excess of those computed from the charac- 
teristics of the landing gear on the assumption 
that the airplane was decelerated a a rigid body 
during the landing impact. It was obvious that 
an adequate stress analysis of the airplane would 
have to take into account the transients excited 
by the landing impact. Several methods were 
soon proposed for computing these transients, 
notably the statistical method of Biot and Bis- 
plinghoff [l], 2 in which the response of the airplane 
for impacts of standard shapes is resolved in terms 
of the normal modes of vibration of the airplane, 
and the method of David Williams [2], in which 



1 Paper presented before Seventh International Congress of Applied 
Mechanics, London, September 1948. 

The work described in this paper was conducted for the Bureau of Aero- 
nautics, Navy Department. 

2 Figures in brackets indicate the literature references at the end of this 
paper. 



the static response and the vibrational response 
are computed separately. A method of coupled 
modes has been suggested recently by Levy, in 
which the modes are computed directly from the 
rigid body motion. Finally, there is a possibility 
of adapting the method of traveling waves that 
was used by St. Venant over 60 years ago to solve 
the problem of longitudinal impact of an elastic 
bar. 

It was decided to try out these methods in the 
computation of flexural transients in an airplane 
wing following a relatively soft and a relatively 
hard, landing impact. It was hoped that the 
computations would indicate the effectiveness of 
these methods for the rapid and accurate compu- 
tation of transients in an airplane wing. 

II. Methods Tried 

1. Method of Normal Modes (Biot and Bisplinghoff) 

The method of normal modes for computing 
the response of an clastic structure to impact is 
described in detail in the classical paper of Biot 
and Bisplinghoff [1], referred to above. The 
method makes use of the independence of the mo- 
tion in one normal mode from that in any other 
normal mode. 
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The total displacement y(x,t) of any point with 
coordinate x is then given as the sum of displace- 
ments VtiXjt) in each one of the modes i that may 
be excited by the impact 



y(x,t)=Jly i (x,t). 

i = l 



(i) 



The displacement yt(x,f) in the ith mode can be 
derived from the displacement q t (t) of the linear 



/////////// 
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Figure 1. Equivalent linear oscillator. 

oscillator of figure 1 in accordance with the equa- 
tion 

y i (x,t)=<p i (x)q i (t), (2) 

where 

<Pi(x)= normalized deflection in ith mode, cor- 
responding to unit deflection of the 
reference station, e. g., the wing tip; 

2<(0= deflection at reference station. 

The deflection q t (t) of the linear oscillator of 
figure 1 under the action of an external force 
Qi(t) is given by the well-known integral of 
Duhamel: 



ii(f)- 



--jhJl Qi{T)sin " i{t ~ T)dT ' (3) 



It is shown in reference [1] that the deflection 
q t (t) has the correct value for the ith mode of the 
elastic structure with its distributed mass m(x) 



under the action of an external impact force 

F(x„t) if 

Mi= f (p 2 i(x)dm= generalized mass in ith^ 
mode ; 



a>i= ■yJK i IM i = frequency of ith mode, 
radians/sec ; 

Qi(t)=<Pi(Xf)F(Xf,t)=geneY&]ized force in 
ith mode. 



(4) 



The stress at any station x is obtained by super- 
posing the stresses corresponding to the deflection 
q t (t) in each mode. 

Numerical solutions for the response q t (t) of a 
linear oscillator to impacts Q t (t) of a variety of 
shapes are given in reference [1]. 

2. Method of Separating Static Deflection and 
Vibrational Deflection (Williams) 

Williams outlined [2] a method of computing 
the displacement of an elastic structure as a 
result of applying an impact force of fixed dis- 
tribution along the structure : 

F(x„t)=P(x)F(t), (6) 

where 

P(z):=spacial distribution of force, assumed 

constant ; 
F\t)= variation of force with time. 
Williams expands P(x) in terms of the inertia 
loads of the various modes 



P(X): 



--i=S.. Lt{x) ' 



(7) 



where L t (x) will be proportional to the product 
of mass and amplitude at each station m(x)<pi(x). 
He reduces the number of natural modes i that 
have to be included in the expansion by noting 
that the inertia forces become negligible compared 
to the elastic force for modes with sufficiently 
high frequency. The deflection in these higher 
modes is therefore the same as if the force were 
acting statically. The deflection in the lower 
frequency modes may be computed by adding a 
static force Li(x)fi(t), equivalent to the inertia 
force, to the external force Li(x)Fi(t). Williams 
shows that the equivalent static force is given by 



»/: 



L i (x)f i (t)=L i (x) cos o)S-T)F(r)dr. (8) 
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The deflection of the structure under the impact 
load is then equal to that resulting from a static 
force, 

P(x)F(t)+XL i (x)J i (f). (9) 

It is interesting to note that a separation of the 
response into a static and a vibrational term can 
be made directly by integrating eq 3 by parts as 
follows : 



*®-ic£«« 



(r) 



COS C0i(t — t)|* 



COf lo 

COS a)i(t — t) 



fW) C0Sa " ( '~ r) ^l 

Jo «< J 



(10) 



Substituting limits and noting that 



gives 
2«(*) a 



''^-J^jyo^S-r)^)^. (11) 



The first term is the deflection of the equivalent 
linear oscillator of figure 1 on the assumption that 
the generalized force Qt(t) acts statically, and 
the second term is the vibrational deflection of the 
oscillator,, 

The resultant response of the structure is given 
by adding to the static deflection under the 
external load F(x f) t) the vibrational deflection 
in each mode. The solution is more general than 
that of Williams in not being restricted to external 
forces with a fixed distribution P{x) in space. 
The two solutions coincide for the special case in 
which F(x f , t) can be written in the form of eq 6. 

3. Method of Coupled Modes (Levy) 

Samuel Levy of this Bureau has suggested that 
it may be possible to represent the response to an 
impact in terms of fewer modes if the restriction 
to normal modes is dropped. He proposed a spe- 
cific series of modes that can be computed directly 
for a structure with known influence coefficients. 

Levy's method will be made clear by applying 
it to the case under discussion of computing the 
flexural transients in an airplane wing following the 
landing impact. 

The rigid-body acceleration y following the 
application of a unit impact force P=l at a dis- 




Figure 2. Force producing rigid-body translation and 
rotation of wing. 



tance x P from the center of the wing in figure 2 is 
given by 

*-i+¥' ' (12 > 

where 

2— distance from center; 
Af =mass of wing; 
/= polar moment of inertia of wing about 
longitudinal axis of airplane. 

The shape <t>i(x) of the first vibrational mode is 
obtained by computing the deflection yi(x) caused 
by the combined action of the applied unit impact 
force and the distributed inertia force yoin(x), using 
the influence coefficients for rigid clamping of the 
wing at the center of gravity. The rigid-body 
motion is eliminated or "swept out" from yi(x), 
since it does not enter in the bending of the wing. 
The elimination is made by subtracting a,i-\-b\X 
from 2/1(2;), where a { and 61 are determined from 
the relations 

[yi(x)—di] m(x)dx=0; 



r. 



f- 



x [yi(x)— bix] m(x)dx-\-0. 



(13) 



<t>i(x) is taken as the deflection resulting after 
normalizing to give unit deflection at the tip on 
the side struck. 

The second vibrational mode </> 2 (x) is derived 
from the first by computing the deflection y 2 (x) 
resulting from inertia forces m(x)<t>i(x) and then 
subtracting a 2 +b 2 x-\-C2<t>i(x) , where a 2 , b 2 and c 2 are 
given by 

I [y 2 (x)—a 2 ]m(x)dx=0, 

x[y 2 (x)—b 2 x]m(x)dx=Q, Y (14) 

fa/2 (x) — #2 — b 2 x — c 2 <pi (x) ] L x (x) dz=0, 
J -1 J 
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and Li(x) is the distribution of load that deflects 
the wing to the shape <t>\{x). <t> 2 (x) is taken as the 
normalized value of y 2 (x) — a 2 — b 2 x— c 2 <t>i(x). The 
last of these equations (14) expresses the condition 
that the load L x {x) should do no work on the 
deflections </> 2 (x). The procedure for the deriva- 
tion of modes of still higher order is obvious. 

The " amplitudes" or generalized coordinates 
qt(t) in the vibrational modes are obtained by 
applying Lagrange's equations. 



d/dT\.dU , ." 



s(I)+ 



dU 

dq 2 



(15) 



where T, U are the kinetic energy and the 
potential energy, respectively, that are stored in 
the wing, and where 4>i(x P )P is the generalized 
force acting on the wing in mode i. The energies 
T, U may be computed by resolving the wing into 
a number of stations. The kinetic energy is then 
given by the sum of the kinetic energies of the 
masses m n at each station n: 



T= 7 z2m n [<pi(x n )q l + <p 2 (z n )q 2 + . . .] 2 . (16a) 



The potential energy is given by the sum of the 
work done by the forces q_\L x { x n ), q 2 L 2 (x n ), ... at 
each station n: 

U=^[qiL l (x n )< Pl (x n ) + qlL 2 (x n )<p 2 {x n )+ . . .] (16b) 

Substituting these expressions for T, U in eq 15 
and carrying out the differentiation will lead to 
equations of the form 

K l q l +M ll q l + M l2 q 2 + . . .=Fi® ' 

K 2 q 2 +M 12 qi + M 22 q 2 + . . .-0, 



in which all quantities are known except the 
generalized coordinates q 1} q 2 , . . . . These equa- 
tions may be solved by writing down separate 
equations for qi(q h q 2 , . . .), §2(21,22 • • .)> etc -> 
and then solving these equations simultaneously 
by a numerical method such as Adams 1 method 
[3, p. 363 to 367]. 



4. Method of Traveling Waves 

It was hoped to derive a method of traveling 
waves for computing flexural transients following 
a sharp impact on an airplane beam. Numerous 
attempts were made to develop such a method, 
but all these failed to lead to a procedure against 
which serious criticisms could not be raised. The 
chief difficulty in this problem as compared to 
St. Venant's problem of the longitudinal waves 
in an elastic bar of constant section seems to 
reside in the basic differential equation of the 
problem. In the case of St. Venant's problem, 
the equation of equilibrium of axial forces leads 
to [4, p. 200, 285] 



where 



dt 2 p dx 2 ~~°> 



(18) 



u (x, t)= displacement along axis of bar; 
E = Young's modulus of material; 
p =mass density of material. 

The general solution of this equation can be 
written in the form 

u=f (V? ~ x ) +fi (Vf +x } (i9) 

where, /, Ji are arbitrary functions, which are 
determined by the initial conditions of the impact. 

The solution given by eq 19 represents two 
traveling waves, one proceeding in the -\-x direc- 
tion and the other in the — x direction. The ve- 
locity of propagation of the waves is -yjE/p, i. e., 
it depends on the properties of the material only. 

A much more complicated differential equation 
applies in the case of flexural waves in an elastic 
bar of constant section. Equilibrium of trans- 
verse forces is expressed in this case, according to 
Timoshenko [4, p. 228], by the equation 



b'y 



&y 



EI^+pA^f-pI 



d.r 



dt 2 



11 Vy 

k'Gbt 



5*y _ EIp 5 4 y 

dx 2 dt 2 

1=0, 



k'G dx 2 dt 2 



(20) 



where 



y(x, £)=lateral deflection; 

EI= flexural rigidity of bar; 
pA=mass per unit length of bar; 

67= shear modulus; 
k f G= ratio of average shear stress to change 
in slope produced by shearing force. 
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The first term in eq 20 is the distributed lateral 
force resisted by simple bending of the beam; the 
seeond is the force required to accelerate the 
sections of the beam in their lateral vibration; 
the third is that required to rotate the sections of 
the beam; and the last two are the forces required 
because of the shearing of the beam. The first 
two terms predominate over the last three in most 
beams: 



«3+* 



dt 2 ' 



<>. 



(21) 



This equation is of a different form than eq 18. 
Equation 20 reduces to a simple wave equation of 
the type of eq 18 only for beams of abnormal 
proportions, in which the second term is negli- 
gible compared to other terms in the equation. 
It follows that a transverse impulse applied to an 
elastic beam of ordinary proportions will not 
travel along the beam without distortion, like a 
longitudinal impulse. 

Nevertheless an attempt was made to solve the 
problem by a numerical method. The beam was 
replaced by a series of concent rated masses con- 
nected elastieally by massless beam elements. It 
was attempted first to find the deflection of the 
masses following the application of a known im- 
pact force to one of the masses by writing eq 21 
as a difference equation and then solving by 
numerical integration for the deflection at each 
station as a function of time. This attempt failed 
because too many stations along the beam were 
required to obtain sufficiently accurate values of 
Eld^y/dx* at each station. 

Next an attempt was made to solve the problem 
on the assumption that the impact was of suffi- 
ciently short duration that the response was 



confined to a small portion of the beam. A small 
number of equations can then be written to express 
the equilibrium between the inertia forces acting 
on the masses in the disturbed portion of the beam 
and the shearing forces transmitted by the beam 
(dements to displace these masses relative to each 
other. This method of attack has been unsuc- 
cessful so far because of the difficulty of keeping 
the extent of the disturbed portion of the beam, 
outside of which all masses can be considered at 
rest, small enough so that only a few influence 
coefficients are required to give the relation 
between shearing forces and displacements. 

III. Solution for Soft Impact and 
Comparison with Experiment 

The first three methods outlined above were used 
to compute the bending strains introduced in a 
model wing following a soft vertical landing impact 
at two points. The impact was regarded as soft 
because the period of impact was about 1.7 times 
the fundamental flexural period of the wing. The 
model wing is shown in figure 3. The wing was 
fabricated from aluminum alloy sheet and angles 
to have a distribution of flexural rigidity and of 
mass (figs. 4 and 5) similar to that of a large four- 
engine transport airplane 4 . A "landing strut," (\ 
was fastened below each one of the two inboard 
engine masses to transmit the landing impact to 
the wing. The landing struts contained a thin- 
walled tube to which wire strain gages were 
fastened to measure the landing impact force as a 
function of time. 

The model was tested by releasing it in a nearly 
strain-free condition (see [5]) from a height of 
about 1 inch to make contact with two synthetic 





Figure 3. Model icing for two-point landing tests 
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"landing fields/' The impact force was recorded 
by the dynamometers in the landing struts, and 
bending strains in the neighborhood of the inboard 
engines were recorded from wire strain gages by 
using the technique described in [5]. The shape 
of the landing-impact force was adjusted by using 
a suitable combination of rubber sheets and mold- 
ing clay for the synthetic landing field. The time 
interval between contact at the two landing struts 
was adjusted by adjusting the relative heights of 
the two landing fields. 

The computations were started by determining 
the first three symmetrical and the first three anti- 
symmetrical flexural modes of the wing. For this 
purpose, the mass of the wing was considered to be 
concentrated at the root and at nine stations along 
each half of the wing. Influence coefficients 
between force and displacement at the mass points 
were computed by treating each half-wing as a 




Figure 4. 



40 60 80 

DISTANCE FROM ROOT.IN. 

Flexural rigidity distribution for model wing 



simple beam clamped at the root. The normal 
modes of the wing were obtained from the influence 
coefficients and from the given mass distribution 
by using a dynamic matrix and iteration procedure 
as explained by Duncan and Collar [6]. The 
deflection at the root of the half-wing (center of 



gravity of the wing) for the symmetrical modes 
was obtained from the condition that the center 
of gravity of the wing must remain at rest for free- 
free vibration. In the case of the antisymmetrical 
modes, the rotation at the root of the half-wing 
was obtained from the condition that the moments 



UJ 

a 



JE 



AAM9 


'III 


r i 


i 


! 


.0002 


- 




- 


.0001 


' i 






- 


O 


— i 1 


.t 


-1 



20 40 60 

DISTANCE FROM ROOTJN. 

(a) 



80 



.03061 



-s. .00673 



.00730 



LB-SEC 2 
IN. 



■O 



J. 



-L. 



-L. 



20 40 60 

DISTANCE FROM ROOT , IN. 
(b) 

Figure 5. Mass distribution for model wing, 

a, Distributed mass; b, concentrated mass. 

about the center of gravity are zero for the free- 
free vibrations. The periods of the natural modes 
are given in table 1 . The first three coupled modes 
according to Levy's method were computed from 
the same influence coefficients as those used in the 
previous computations, following the procedure 
described in section II, 3. 

Table 1. Periods of computed natural flexural modes 



Mode 


Period 


First symmetrical ... 


Seconds 

0.0315 

.0115 

.00613 

.0170 
.00759 
.00500 


Second symmetrical . 


Third symmetrical 


First antisymmetrical. 

Second antisymmetrical 

Third antisymmetrical 
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The response of the wing in each mode accord- 
ing to the normal-mode method and according to 
Williams' method was then computed by substi- 
tuting the normalized modes and the generalized 
force obtained from the measured impact force at 
each landing strut in eq 2, 3, 8, and 9. The 
response according to Levy's method was obtained 
by solving eq 17 as indicated in section II, 3, sub- 
stituting the given impact forces at the two land- 
ing struts and then superimposing the two solu- 
tions. 

The bending moments in a wing section near 
the two landing struts were computed as a func- 
tion of time from the amplitudes q t in each one of 
the component modes. 

The resultant bending moments thus obtained 
are plotted against time in figure 6. The figure 



also shows the bending moment corresponding to 
rigid body motion and the observed bending 
moment obtained from measurements of extreme 
fiber-bending strain in the two-point landing test. 
It is apparent that the maximum values of the 
bending moments due to rigid-body motion 
amounted to about 85 percent of the maximum 
values of experimentally measured bending mo- 
ments. The measured bending moments were in 
close agreement with those obtained by Levy's 
method and by Williams' method. It can be seen 
from figure 6 that the difference in maximum 
measured bending moments and maximum com- 
puted bending moments was less than 5 percent 
for Levy's method and less than 7 percent for 
Williams' method. For the method of normal 
modes, this difference between experimental and 
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Figure 6. Response for soft, unsymmetrical, two-point landing impact. 

a. Beading moment, 1 in. outboard of landing strut hitting first; b, bending moment, 1 in. outboard of landing strut hitting second. 
O, Observed; N, normal modes method; W, Williams' method; L, Levy's method; R, rigid body method. 
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Figure 7. Components of response for soft, unsi/m- 
metical, two-point landing impact, normal modes method. 

a, Bending moment, 1 in. outboard of landing strut hitting first; b, bending 
moment, 1 in. outboard oflanding strut hitting second. 

N, Normal modes method (sum of the following components); si, first 
symmetrical mode; s2, second symmetrical mode; s3, third symmetrical 
mode; al, first antisymmetrical mode; a2, second antisymmetrical mode; a3, 
third antisymmetrical mode. 

computed maxima was as much as 23 percent. 

The reason for obtaining better agreement by 
Williams' method than by the method of normal 
modes becomes apparent from an examination of 
figures 7 and 8, which show the contributions to 
the bending moment of the six different modes 
used in the computations. The normal modes 
method, figure 7, shows good convergence for the 
symmetrical modes; however, the third antisym- 
metrical mode still contributes, 13 percent to the 
measured maximum bending moment. Higher 
modes should probably be considered. In the case 
of Williams' method, figure 8, the third antisym- 
metrical mode contributes only about 1 percent 
to the measured maximum bending moment. 
The use of the first two symmetrical and the first 
two antisymmetrical modes in Williams' method 
would have sufficed to give the maximum bend- 
ing moment within 8 percent. 

It is difficult to judge convergence of Levy's 
method by comparing the response in the higher 
modes to the total response, since in this case the 
modes are coupled and the addition of a higher- 
frequency mode term changes the response in the 
lower-frequency modes. Hence the change in 



computed maximum bending moment due to the 
addition of a higher-frequency mode may be ap- 
preciably less than the response in that mode. 
Nevertheless, it is of interest to compare the 
response in the three modes. Figure 9 shows this 
response. The third mode contributed about 15 
percent to the measured maximum bending 
moment. 
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Figure 8. Components of response for soft, unsyinmetrical, 
two-point landing impact, Williams' method. 

a, Bending moment, 1 in. outboard oflanding strut hitting first; b, bending 
moment, 1 in. outboard of landing strut hitting second. 

W, Williams' method (sum of the following components); R, rigid body 
method; si, first symmetrical mode; s2, second symmetrical mode; s3, third 
symmetrical mode; al, first antisymmetrical mode; a2, second antisym 
metrical mode; a3, third antisymmetrical mode. 

IV. Solution for Hard Impact 

The hard impact was chosen as identical in 
shape to the soft impact, but the duration of the 
impact was reduced to one-fifth of that of the soft 
impact. The period of the impact was about 0.34 
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times the fundamental flexural period of the wing. 
The response to the hard impact was computed 
keeping all quantities identical with those for the 
soft impact, except, of course, the duration of 
impact. 

The resultant bending moment near the inboard 
engine above the landing strut, making first con- 
tact with the landing field, is given in figure 10. 
The figure also shows the bending moment corres- 
sponding to the rigid body motion of the wing. 

It is apparent that there is a larger divergence 
between the results of the three methods than for 
the soft impact and also that there is a considerable 
lag between the resultant response and the rigid 
body curve. The maximum bending moment 
according to the method of normal modes is 35 
percent below that for the rigid body motion; 
that for Williams' method is 3 percent below, and 
that for Levy's method is 19 percent above. 

The response in the individual modes of vibra- 
tion according to the three methods is given in 
figures 11 to 13. The convergence is good only for 
the symmetrical modes. The third antisymmetri- 
cal mode according to the method of normal modes 
contributes as much as 23 percent to the computed 
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Figube 9. Components of response for soft, unsymmetrical, 
two-point landing impact, Levy's method. 

a, Bending moment, 1 in. outboard of landing strut hitting first; b, bending 
moment, 1 in. outboard of landing strut hitting second. 

L, Levy's method (sum of the following components); 1, first mode; 2, 
second mode; 3, third mode. 
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Figure 10. Response to hard, unsymmetrical, two-point 
impact. 

Bending moment, 1 in. outboard of landing strut hitting first. 
N, Normal modes method; W, Williams' method; L, Levy's method 
R, rigid body method. 

maximum bending moment (fig. 11). For Will- 
iams' method, the third antisymmetrical mode 
contributes 15 percent to the computed maximum 
bending moment (fig. 12). For Levy's method, 
figure 13, the third coupled mode contributes 43 
percent to the computed maximum total bending 
moment. It is apparent from this that a fourth 
analytical procedure needs to be developed to de- 
rive the response to impacts of very short duration, 
such as that chosen in the example. 

V. Conclusions 

Analysis of the flexural transients in a model 
airplane wing following a soft unsymmetrical two- 
point landing impact indicates that Williams' 
method, separating static and vibrational response, 
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Figure 11. Components of response for hard, nonsym- 
metrical, two-point impact, normal modes method. 

Bending moment, 1 in. outboard of landing strut hitting first. 

N, Normal modes method (sum of the following components); si, first 
symmetrical mode; s2, second symmetrical mode; s3, third symmetrical mode; 
al, first antisymmetrical mode; a2, second antisymmetrical mode; a3, third 
antisymmetrical mode. 



-200 



.004 .008 .012 .016 

TIME .SECOND 

Figure 12. Components of response for hard, unsym- 
metrical, two-point impact, Williams 1 method. 
Bending moment, 1 in. outboard of landing strut hitting first. 

W, Williams' method (sum of the following components) ; R, rigid body 
method; si, first symmetrical mode; s2, second symmetrical mode; s3, third 
symmetrical mode; al, first antisymmetrical mode; a2, second antisym- 
metrical mode; a3, third antisymmetrical mode. 
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Figure 13. Components of response for hard, unsy in- 
metrical, two-point impact, Levy's method. 

Bending moment, 1 in. outboard of landing strut hitting first. 
X, Levy's method (sum of the following components); 1, first mode; 2, 
second mode; 3, third mode. 

is superior to the normal modes method used by 
Biot and Bisplinghoff , in leading to values in good 
agreement with experiment without using more 



than four modes. A method of coupled modes 
suggested by Levy leads to values of the response 
about as accurate as those from Williams' method 
with as few as three vibrational modes. 

None of the three methods showed rapid con- 
vergence for the response following a hard impact 
of the same shape as the soft impact, but with only 
one-fifth of the duration. Attempts to develop a 
traveling-wa^e method for the solution of this 
problem have been unsuccessful so far. 



The author expresses his indebtedness to his 
colleagues Samuel Levy, J. B. Woodson, and 
W. D. Krol], who were principally responsible for 
the computation of the transients given in the 
paper. 
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